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ROLE OF GAS-SURFACE INTERACTIONS IN THE REDUCTION OF 
OGO 6 NEUTRAL-PARTICLE MASS SPECTROMETER DATA 


A. E. Hedin 

Goddard Space Flight Center 

B. B. Hinton and G. A. Schmitt 
University of Michigan 


INTRODUCTION 

Several authors have pointed out that gas-surface interactions complicate the interpretation of 
neutral-particle density measurements in the thermosphere (Moe and Moe, 1967, 1969; Silverman and 
Newton, 1970; von Zahn, 1967). This paper describes the gas-surface interaction effects observed by 
a satellite neutral-particle mass spectrometer and summarizes the technique developed to account for 
them in determining ambient neutral densities. 

DESCRIPTION OF THE EXPERIMENT 

The Orbiting Geophysical Observatory (OGO) 6 is a polar-orbiting satellite with an inclination of 
82° prograde, an apogee of approximately 1100 km, and a perigee of about 400 km. It was launched 
on June 5, 1969 (day 156 of 1969), and was operated until July 1971. The neutral-particle mass 
spectrometer is mounted in an orbital-plane experiment package (OPEP) such that in the normal mode 
of operation, the spectrometer inlet faces the spacecraft’s direction of travel. In this mode, the flux of 
ambient neutral molecules entering the spectrometer is modulated primarily by the varying ambient 
density as the spacecraft changes altitude and position. On command, the OPEP can be rotated 
cyclically through 220° about a vertical axis, effecting an angle-of-attack modulation of the incoming 
flux of ambient molecules. 

The quadrupole mass spectrometer employs an electron-impact ion source behind a gold-plated 
stainless-steel antechamber (fig. 1 ). Ambient neutral molecules enter the antechamber through a knife- 
edged orifice and undergo many surface collisions before entering the ion source region. Ions created 
by impact with a beam of 75-V electrons in the source region are focused by an electrostatic lens and 
filtered according to their charge-to-mass ratios. Details of the instrument construction and operation 
have been documented elsewhere (Carignan and Pinkus, 1968). 

In the normal mode, the total ion current and the ion currents for ions with molecular weights 2, 
4, 16, 28, and 32 are sampled for 1.152 s once every 9.216 s during 258 s out of a 368-s cycle. During 
1 10 s of the 368-s cycle the spectrometer is in sweep mode, during which period each ion in the range 
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Figure 1. -Cross section of mass spectrometer antechamber 
and ion source. 


from 1 to 46 amu is sampled at least twice. A continuous-sweep mode is available on command, and 
masses of 2, 4, 28, 32, or the total flux can be sampled continuously on command. 


GENERAL OBSERVATIONS 


Based on the assumption that incoming ambient molecules are thermally accommodated to the 
antechamber surface temperature, the following equation gives the number density as a function of 
time t for any constituent in the ion source region: 


dn g _ (1/4 )n a F{s)c a A h - (l/4)n g cA h 
dt V 

where 

2ir 1 l 2 v n 

F(s) - for v n ^ > c a (2) 
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and 


v n = spacecraft velocity component normal to the inlet orifice 
n - ion source number density 
n -ambient number density 
A h = area of the inlet knife-edged orifice 

c = mean thermal speed of constituent molecules at the antechamber surface temperature 


c a = mean thermal speed of constituent molecules in the ambient gas 
V - volume enclosed by the antechamber 

Equation (1) describes the balance of molecular fluxes into and out of the spectrometer and may 
be used to consider the response of the source density to changes in incoming flux. For example, if 
the incoming flux decreases to zero, then 


tog = c *h 
dt 4V 


(3) 


and the source density would decrease exponentially with a time constant of 4V/cA /j . This value for 
the OGO 6 mass spectrometer is of the order of 0.02 s. Physically, this means that the antechamber 
should be emptied of molecules almost immediately if the incoming flux is stopped. 

A direct test of the simple model represented by equation ( 1 ) is provided by a comparison of 
figures 2(a) and (b). Figure 2(a) is a plot showing typical source-density data taken during one orbit 
of the OGO 6 spacecraft. The observed source densities plotted in figure 2(a) were obtained by 
multiplying the ion current for each constituent by the laboratory sensitivity for that constituent. 
Figure 2(b) is an equivalent plot using equation (1) together with a model atmosphere prediction of 
ambient densities as a function of altitude. 

Figure 2(a) indicates that atomic oxygen has combined in large part to form molecular oxygen, 
which appears as an ion current at mass 32 and a smaller ion current at mass 16. Mass 16 values are 
largely the result of fractionation of molecular oxygen in the source region. Both mass 32 and mass 28 
have a maximum source density near perigee as expected from the strong dependence of neutral 
ambient density on altitude. There is, however, an asymmetry in the source density about perigee due 
to surface absorption processes before perigee and desorption processes after perigee. All the con- 
stituents except helium have detectable source densities near apogee as a result of desorption. Finally, 
the presence of mass 44 (C0 2 or N 2 O) and mass 1 (atomic hydrogen), as well as carbon monoxide, is 
thought to result from chemical reactions of atmospheric atomic oxygen with constituents on the 
instrument surfaces. 

Figure 3 shows some long-term features of the measured source-density data. The perigee and 
apogee source densities during an orbit are plotted versus time for ions with masses 16, 28, and 32. 

The perigee density data reflect the changing atmosphere as the spacecraft perigee changes latitude and 
local time, and as the atmosphere itself varies. The apogee source density is higher than simple theory 
predicts (fig. 2(b)) and depends to some extent on the magnitude of perigee density. It was also found 
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Figure 2.-Source densities during orbit, (a) Observed, for typical orbit, (b) Hypothetical, without gas-surface 

interactions. 


that, after an initial period, the maximum densities for masses 16 and 32 have a nearly constant ratio 
of 0.1, indicating a common source for these two constituents and a certain stability in the processes 
leading to these signals. Also, the long-term trend in minimum density for mass 28 is decreasing, 
unlike mass 32, indicating that the processes determining minimum densities are not the same for these 
two components. 

Several problems arose during the operation of the experiment that affect understanding of the 
gas-surface interactions in the spectrometer. At approximately 30 s after turn-on there was a sudden 
loss of sensitivity by approximately a factor of 5 (5.03 for mass 32, 4.70 for mass 28, 1 .82 for mass 
16, 4.32 for mass 4, and 3.46 for mass 2). The exact cause was never determined but was believed to 
have originated in arcing nearby or within the experiment electronics. No indication of a further 
change in sensitivity was found and all densities have been adjusted for this sensitivity change. The 
change in sensitivity was nearly the same for every mass except 16, and the mass 16 to mass 32 source- 
density ratio before the shift was near the value expected from fractionation (0.0455). Thus it must 
be assumed that the fractionation ratio of molecular oxygen increased at this instant to 0.16. Later, 
however, the 16 to 32 ratio was always lower than this value. Despite various laboratory tests on a 
prototype instrument, this failure pattern could not be reproduced and, therefore, the detailed 
mechanism is not clear. The observation of a constant 16 to 32 ratio after the first 2 weeks indicates 
that the new sensitivities have been stable. Nevertheless, the difficulty in determining the correct value 
to be used for the fractionation pattern of oxygen contributes to the inability to uniquely determine 
certain surface parameters. 
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On day 225 of 1969, an electronic failure occurred which put the quadrupole spectrometer into 
an anomalous operating mode during the measurement of mass 32. The net result was that the ion 
current recorded for mass 32 became the total ion current for masses greater than 27. To compensate 
the data acquired after this failure, the mass 28 ion current was subtracted from the mass 32 ion 
current in the analysis process. Unfortunately, there were some small currents at other masses, partic- 
ularly 30 and 44, which could not be determined with the frequency or accuracy needed to permit 
subtraction from the measured mass 32 ion current. The mass 32 source density measurement thus 
contained a residual contribution from these other masses, as seen by the jump in apogee value at day 
225 in figure 3. The principal effect on ambient density determination is not from the direct con- 
tribution, which is negligible near perigee, but in the error it may cause in surface parameter 
determination. 

The temperature of the antechamber is of interest because the surface desorption rate is an 
Exponential function of temperature (Ehrlich, 1959). A thermistor, located near the source region 
(fig. 1), recorded changes in temperature of ±10 K with a nominal value of 293 K. These temperature 
variations are sufficiently large to produce noticeable changes in desorption— nearly a factor of 4 for a 
desorption energy of 92 kJ (22 kcal) per mole. However, the measured temperature may not 
accurately reflect the effective antechamber temperature because, for example, the measured tempera- 
ture was minimized when the angle between the sun and the front of the antechamber was smallest- 
opposite to the anticipated behavior. The measured temperature was ignored in the final surface 
analysis, but these data were examined for possible examples of thermal effects. 

MOLECULAR NITROGEN DENSITY 

Analysis of gas-surface interactions affecting the mass 28 source density depends critically upon 
whether the ion current at mass 28 is due primarily to N 2 or CO. If the mass 28 ion current is entirely 
due to N 2 , then one could use an adsorption model for N 2 of the type discussed by Moe and Moe 
(1967, 1969) in connection with the analysis of pressure-gage data. The net result will be an additive 
correction before perigee and a subtractive correction after perigee. If, on the other hand, the mass 
28 ion current at apogee is due primarily to CO, one must consider the adsorption-desorption 
reactions of this gas rather than N 2 , and the correction will always be subtractive. 

The most direct determination of apogee mass 28 peak composition would come, in principle, 
from comparing the mass 14 and 28 ion currents. Mass 14 ion current results from dissociative 
ionization of N 2 , while mass 28 ion current results from both N 2 and CO. Early in the spacecraft’s 
lifetime, the ion current measured at mass 14 indicated that CO was present near perigee. During this 
early period there were substantial ion currents at mass 44, indicating the probable presence of neutral 
C0 2 . Unfortunately, at high altitudes where the concentration of N 2 in the atmosphere is negligible, 
the ion current at mass 14 was too small to identify the neutral constituent contributing to the mass 28 
ion current. The drop in sensitivity referred to earlier resulted in a generally undetectable mass 14 ion 
current at apogee after the first few days. If the apogee mass 28 ion current were entirely due to N 2 , 
the expected mass 14 ion current would have produced only one or two bits of telemetry signal. 

However, the identification of the mass 28 background gas as CO is supported by other evidence. 
Laboratory experience indicates that CO is almost always present to some degree in evacuated systems. 
The long-term downward trend in the apogee mass 28 source density is not compatible with the apogee 
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density resulting from N 2 desorption, unless an ad hoc change in the desorption rate coefficient is 
postulated. Furthermore, there was a definite, measurable increase in CO at filament change on day 
286 of 1969, thus indicating the continuing importance of CO. Finally, previous laboratory work on 
surface interactions (Trapnell, 1953) indicates that gold is the least reactive of all metals and N 2 the 
least reactive molecule (except for the noble gases) on any metal, whereas CO is known to chemisorb 
on gold. 

To develop a definite model, one may suppose that there exists an initial surface density of 
carbon that is slowly replaced by diffusion as it is depleted by reaction with oxygen atoms. The base 
metal, stainless steel, would be a reservoir for carbon that could slowly diffuse through the base metal 
and the gold-over-nickel surface layer. Extrapolation of available data indicates that the total CO 
generated in the instrument will consume a number of carbon atoms equivalent to at least a few parts 
per million of the antechamber base metal. 

It is assumed that CO is produced by two reactions involving wall w and gas g atoms: 

C(w) + 0(g) ^CO(g) 

C(w) + O(w) CO(vv) 

Further, it is assumed that CO is destroyed by three reactions: 

C0(w) + 0(g)->C0 2 (g) 

C0(g) + 0(w)-+C0 2 Cg) 

CO(w) + O(vv) -» C0 2 (g) 

Although these reactions represent negligible sinks for oxygen atoms (after the first 20 days of flight), 
they do constitute a significant loss of CO. 

In writing these equations, it is assumed that some of the incoming oxygen atoms are adsorbed 
on the walls. Molecules (or atoms) on the wall are distinguished as a separate species (e.g., CO(w)) 
from the ordinary, or gas phase molecules or atoms (e.g., CO(g)). It has been assumed that there are 
at most two adsorbed states of CO and a single adsorbed state of atomic oxygen. It is thought that 
most chemisorption systems have multiple binding states, each with a distinct sticking coefficient (or 
probability) and desorption time constant. It is conceivable, for example, that the CO molecules 
adsorbed from the gas phase are bound to the surface mainly in one state and the molecules produced 
by adatom-adatom combination in a preadsorbed state might be bound primarily in a second state. 

Any attempt to account for this and for the rates of interconversion among the states increases 
enormously the complexity of the problem. These considerations would be appropriate for laboratory 
investigations where, in principle at least, conditions may be controlled and are repeatable. As far as 
possible, all but the dominant states have been ignored because of the lack of sufficiently detailed 
information to consider more than one or two surface states. 

The following expressions were obtained by considering production, loss, outflow, adsorption, 
and desorption. In these equations and in those in the following sections, k is used to denote sticking 
coefficients, a to represent cross sections for reactions of an adsorbed constituent with one in the 


7 



gas phase, o aa to represent a cross section per unit time that describes the rate of reaction of two 
adsorbed constituents, and t w to represent the lifetime of an adsorbed constituent. A subscript w on 
a density refers to the number per unit area on the wall; a subscript g refers to an ordinary volume 
number density. All of these quantities are associated with a specific process by means of additional 
subscripts and superscripts, as indicated in the following list: 

A h = inlet aperture area, cm 2 

A w - wall area of antechamber, cm 2 

cj 6 = mean thermal speed of ambient atomic oxygen, cm s _1 
cj 6 = mean thermal speed of atomic oxygen, cm s -1 

c 28 = mean thermal speed of molecular nitrogen or carbon monoxide, cm s" 1 
cf = mean thermal speed of molecular oxygen, cm s -1 
& 16 = sticking probability of atomic oxygen 
k CQ - sticking probability of carbon monoxide 

m] 6 = density of ambient atomic oxygen multiplied by the function F(s), cm -3 
n j 6 = density of atomic oxygen in antechamber, cm -3 
2 = density of molecular oxygen in antechamber, cm -3 
n g° = density of carbon monoxide in antechamber, cm -3 
n ^ 2 = density of carbon atoms on walls, cm -2 
n = density of atomic oxygen on walls, cm -2 
= density of carbon monoxide on walls, cm -2 
V= antechamber volume, cm 3 
°aa = a datom-adatom recombination coefficient, cm 2 s _1 
°ag = adatom-gas atom recombination coefficient, cm 2 
a aa° = coefficient of CO production by C(w) + O(w), cm 2 s -1 
Og® = coefficient of CO production by C(w) + 0(g), cm 2 
°a° 2 = coefficient of C0 2 production by CO(w) + 0(w), cm 2 s _1 
a< ag 2 ~ coefficient of C0 2 production by CO(w) + 0(g), cm 2 
°ga° 2 ~ coefficient of C0 2 production by CO (g) + 0(w), cm 2 
1 /r^ 6 = desorption rate of atomic oxygen, s -1 
1 /tJ = desorption rate of carbon monoxide, s _1 
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Outflow is sufficiently rapid that the steady-state assumption applies to n^° . Consequently 


dn co nl 6 cl 6 n} 2 A, 


n c0 A, 


n co c 28 


1- = " S V co + ( , co/4w + i4fc + ^o 2n i6 Aw)=0 (4) 


a w + 

- Fr co 


dn c .° 


k~~rr; u c 


C0^28 


w i t rn i & & CO 

— =n 1 ?(£?n}f + - 


w aa w 


w \ T co 


+ a co 2« 16 + 
aa w 


a co 2 c 16 n 16 
ag g g 


Rearrangement of equation (4) results in 


4 „CO /t CO + CO 16^16 12 

„co ^ Iz <* » ± " - 4 (6) 

‘ 2 “(*co K+A„+oC°inlU w ) " 

if dn co Idt is negligible in comparison to the other terms in the equation. If one considers n l g 6 and 
n to be explicitly known functions of time, equation (5) can be solved for n^° if n g is replaced by 
equation (6). If one lets 

1 i , ,, ^co A w 

r (*\ = j- fyC0 2 + — qC 0 2 ^,16^16 _ " (7) 

t co " “■ 4 « * ' rl°ik co A w *A h *of°2n'J'A w ) 


«o=< : k>l 6 + 


„16t: 16 a „CO 
CO U g C g A »°ag 

4 k C0 A w + A h+a^ 0 2nl 6 A Wi 


+ n<=° L{t) = Pit ) 


and the solution is 


° = exp 


- r ) dt' f 

J _ oo •/ — oo 


P(f) exp I 1(f ) dt dt + c exp 




where c is a parameter, discussed later. After sufficient time 


«c° = exp 


T-f £(*')#' T 

I J — oo %/ -oo 


/»(/') exp / L(t")dt" dt' 


Use of this result with equation (6) yields 


o^ 0 ril 6 cl 6fl l 2 A w + (4A W /r£° ) ]p * exp / L(r')df 


_C0^28 /j. a + a +a co 2n 16 A ) 
r w ^ 71 ft V n w^w> 


in which * denotes convolution. 
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In the investigation of CO, using a heuristic approach, it was discovered that the CO density could 
be approximated by an expression of the form 


,co 



32 

n g exp 



dt' + c 


(13) 


The CO density during an orbit was determined by subtracting the mass 14 source density from the 
mass 28 source density. The mass 14 source density took into account the fractionation ratio and 
should have agreed with the mass 28 source density if no CO were present. Also, the mass 28 source 
density within 1500 s of apogee was taken to indicate CO density, consistent with the earlier decision 
on the composition of the apogee mass 28 density. After several months, the mass 14 source density 
near perigee was consistently higher than the mass 28 source density by about 1 1 .4 percent. There 
was no way to determine from the available measurements whether this was due to an error in the 
fractionation ratio or to a true, ambient, atomic nitrogen density. It was decided to adjust the frac- 
tionation ratio to make the two source densities agree during the later orbits. 


Equation (13) can be interpreted in terms of a simplified solution of equation (12) in which 
a oa° > °aa° 2 > a ag° 2 > an< * °g a ° 2 are ignored and n 22 is assumed to be a rough measure of nj 6 . The result- 
ing solution for n^° is given by equation (14), in which k’ QQ = k CQ A w /A h , a modified sticking 
coefficient, r£° = (1 + k' CQ )r£° , a modified desorption time constant, and a = a' (l + k ' CQ ), a param- 
eter proportional to n^ 2 : 


„CO 


,32 


Ir' 

K co 


1 +k' 


CO 


1 + k' 


CO 



+ c 


(14) 


Only during the first 2 weeks of flight were there sufficient CO data near perigee to allow a 
determination of all four unknown parameters. The data were grouped into half-day intervals, and 
parameters a, , and c were determined for various fixed values of r^ 0 ’ . A value of 1 500 s for r co ’ 
provided the best fit to the data. An example of the type of fit obtained is shown in figure 4, 



Figure 4.— Measured CO during a typical early orbit 
compared with model estimates. 
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Figure 5.-Variation of parameter A^q for the period when CO density 
was available near perigee. 


including the effect of measured antechamber temperatures. Because temperature variations over the 
early orbits were small, it made little difference whether they were included or not. After launch on 
day 1 56, the value of k' CQ varied daily as shown in figure 5. It appeared to level off to a value of 
about 2.5 by day 180 and was impossible to determine after that. After day 180, both t£° and k QO 
were held fixed; parameters a and c were allowed to vary, and values were then determined from all 
the data in a day. Thus, after day 180 it was assumed that 


,co 


r »i 2 


2.5 


L3.5 3.5 



(15) 


It is probable that the parameter c in equations (13) through (15) represents a process with a time 
constant long compared to the orbital period. That is, c is not a true constant, but might represent a 
quantity that changes very little around an orbit. This leads one to suppose that two surface states 
might be significant for CO, with c representing the more slowly changing one and the expression 
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the faster one. If this is the case, each surface state must be considered as a separate species governed 
by its own set of equations having the form of equations (5), (9), (10), and (11). After the appro- 
priate modifications to equations (4) and (6), the solution for of equation (12) is replaced by 


c ] g 6 ng 6 n l w 2 A h o%p + 4A W { P * exp 


»co 


-f L(t')dt' 

J — 00 


/r £0 + T> * exp 


f L\t')dt' /, 

_ J — oo _ 


! r co # 

* w 


(^“CO + + ^ h + °ag 


(16) 
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The symbols P' and L' refer to the second state and are defined by relations analogous to equations 
(7) and (8). 

Equation (15) may be regarded as an approximate expansion of equation (16) in terms of the 
convenient function n 22 . To the extent that the shape of n g 6 resembles n 22 , one may identify the 
first term of equation (15) with the first term of equation (16), the second term with the second, and 
so on. 

The long-term variation of parameter a (proportional to the surface wall density of carbon atoms) 
can be predicted under the assumption that this density is diffusion controlled. If n 12 arose from a 

W 

diffusion of carbon atoms through the wall of the antechamber, then one would expect that n\ 2 would 
vary as exp (~ir 2 Dt/d 2 ), where D is the diffusion coefficient and d is the wall thickness. The diffusion 
coefficient is not known precisely, but is of the order of 3 X 10 -10 cm 2 s _1 at 300 K (Kittel, 1968). 

Figure 6 is a plot of the a values for the first 6 months of spacecraft operation. The curves marked 
A represent empirical Fits to the data points of the form 

/-Dy\ /-Dy\ 

a = a i exp 1—1 +a 2 exp I — ) (17) 

where a l ,a 2 ,T l , and T 2 are constants and Dy is the day number. Because of the obvious perturba- 
tions occurring at the filament change on day 286, a second curve was used beginning on that day. 

The B curve has been calculated for comparison, using a diffusion coefficient of 3 X 1CT 10 cm 2 s -1 , 
accounting approximately for the various wall thicknesses of chamber parts. 

The most obvious feature of the plot shown in the figure is the cyclic variation of parameter a. 

This is not consistent with the conception that a is proportional to the carbon available for reaction 
with the incoming oxygen. One expects this parameter to decrease monotonically with time and its 
profile to have a greater slope when ambient oxygen is highest. Thus, either there must be some ex- 
planation for a carbon variation of this form, or the apogee slope variation must have a contribution 
from some source not included in the model. The large excursions of a also reflect data-analysis prob- 
lems, because a was determined primarily from apogee data with a low signal-to-noise ratio. Also, the 


Figure 6.— Variation of parameter a found from apogee 
data. Curves A are empirical fits to data and reflect the 
ion source filament change. Curve B results when carbon 
surface-density is limited by diffusion through the thin 
gold film. 
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(c) TIME FROM PERIGEE ( s) 



Figure 7-Mass 28 and CO source densities, 
showing temperature effect, (a) Data near apogee, 
(b) Near-apogee data for different orbit, (c) Data 
for an entire orbit. 


excursions might be artifacts of the approximations leading to equation (16). And, finally, the phe- 
nomena in question are extremely temperature dependent. 

The effect of variations in surface desorption rate resulting from variations in surface temperature 
was studied by recalculating the a parameter, using temperature as measured by the source thermistor. 
There was no basic change in the long-term variation of a. Also, individual orbits were examined to 
clarify the effect of including temperature. A plot of mass 28 source density near apogee is shown in 
figure 7(a) compared with the best-fit prediction determined both with and without the temperature 
variation of desorption. Note that the sign of dnjdt can be reversed by including temperature varia- 
tions. While the data have considerable scatter, the average behavior is better fitted by including the 
temperature variation. On the other hand, figure 7(b) shows an orbit where inclusion of the tempera- 
ture variation makes the fit worse. Figure 7(c) shows the difference over the whole orbit and indicates 
that there is little practical difference near perigee. The ambiguity of these results may be due to the 
placement of the thermistor on the outer housing of the antechamber-ion-source assembly. Although 
this gives a general indication of the temperature in the relevant region, it is probable that temperature 
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differences of the order of 2 or 3 occur, and that the phase of temperature variation around the 
orbit at the thermistor might be somewhat different than that inside the antechamber itself. 

There is a suggestive correlation between the period of the variation in the a parameter and the 
variations in many spacecraft parameters as a result of the slow (2° /day) drift of the orbit plane with 
respect to the Earth-Sun line. If one considers a given point in the orbit, say apogee, the angles be- 
tween the spectrometer orifice and the Sun or solar panels and the measured source temperature all 
have cyclic variations similar to that of the a parameter. 

A possible cause of variation in density slope near apogee would be detection by the spectrometer 
of outgassing from the solar panels. Two arguments against this are the lack of detectable variation in 
the mass 28 peak at apogee during OPEP scans and the fact that the minimum value of parameter a 
occurs when the orifice looks at the body before apogee and away from the body after apogee. The 
latter condition is always associated with noon-midnight orbits. 

Figure 8 shows a typical orbit, comparing the observed mass 28 source density with model esti- 
mates of the contribution of CO to that source density. Figure 9 shows the fractional error in calcu- 
lated molecular nitrogen ambient density due to errors in a for the data in figure 8. Errors are com- 
pared with results from full surface analysis. Note that errors in a affect calculated ambient density 
more severely after perigee and that the calculated ambient density is relatively insensitive to large 
excursions in a. This effect is the inverse of the sensitivity of a to errors in the data. 

It is clear from figure 4 that the CO model described is capable of giving a reasonable fit to the 
measured CO density over an orbit. The problem arises in later orbits when CO density near perigee 
becomes small, with large errors, and only the apogee density can be used in determining parameters 




Figure 8.-Comparison of mass 28 source density with predicted Figure 9.-Sensitivity of calculated N 2 ambi- 

CO density for several values of a. The calculated variation of ent density to changes in surface parameters, 

atomic oxygen wall density and measured source densities for 
masses 32 and 28 are shown. 
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of the model. The small CO density at perigee means that the correction is likewise small at perigee, 
but the CO prediction becomes critical at higher altitudes and its accuracy is one of the factors that 
limits the altitude of reliable N 2 densities. Because of the danger that variations in the a parameter may 
be the result of extraneous effects not included in the model, a smooth curve was drawn through the 
values, as indicated in figure 6, and a 50-percent error bar was imposed on the CO correction (inter- 
preted as the difference between the predicted CO and the apogee CO) after day 180 of 1969. For the 
earlier orbits, the CO correction was taken to be accurate to 25 percent. For reliable estimates of N 2 
densities, it appears that a reasonable correction for CO can be made with the method described and 
that the use of N 2 density estimates with large corrections can be avoided. 

ATOMIC OXYGEN DENSITY 

A model of the atomic oxygen interactions in the OGO 6 mass spectrometer has been constructed 
that accounts for the general behavior of the data. The model is a quasi-equilibrium model that ignores 
density gradients in the antechamber and variations in surface properties. It is applicable only after 
approximately 20 days of flight when reactions producing CO and C0 2 become an insignificant loss 
process for oxygen. The model assumes that the flux of particles into the antechamber from the 
ambient atmosphere is primarily atomic oxygen with an insignificant component of molecular oxygen. 
The principal processes are expected to be adsorption of atomic oxygen and gas atom-adatom recom- 
bination, with a small amount of adatom-adatom recombination (Wood, Baker, and Wise, 1970). 
Desorption of atomic oxygen is included to explain the finite atomic oxygen source density at apogee. 
The interaction of molecular oxygen with the surfaces is ignored, and recombined oxygen is assumed to 
appear immediately in the gas phase, because molecular oxygen is not expected to chemisorb on gold 
(Trapnell, 1953; Wood, Baker, and Wise, 1970) and would thus have a very short surface interaction 
time. 

The rate of adsorption of atomic oxygen on the antechamber surfaces is proportional to the rate 
at which gas atoms strike the surface. From elementary kinetic theory, the collision frequency per 
unit area of a gas with volume number density n and mean thermal speed c is n g c/ 4. If n represents 
the number of potential sites (usually taken to be 10 15 cm -2 ) and n w the number of sites that are 
occupied, then the proportionality constant is (1 - n w /n)k, where k, the sticking coefficient, is 
regarded as constant (Langmuir isotherm). It was found that the oxygen surface coverage, using 
parameters that gave reasonable agreement with the data, was always of the order of 10 14 cm 2 or 
less, and thus the factor (1 - n w /n ) was taken equal to 1. 

The rate of gas atom-adatom recombination is proportional to the collision rate of gas atoms with 
the wall and to the number of adsorbed atoms per unit area. The rate of formation of molecules is 
thus o a n n w c/4, where o ag is the recombination coefficient. The rate of adatom-adatom recombina- 
tion is proportional to the square of the wall coverage, and thus the rate of formation of molecules is 
given by a (n w ) 2 , where o aa is the recombination coefficient. The adatom-adatom recombination 
rate is a function of temperature such that o aa is proportional to exp ( -E/RT ), where E is the energy 
required for surface diffusion and reaction activation, R is the gas constant, and T is the absolute 
temperature. The desorption of surface atoms occurs at a rate given by n w /r, where r is a desorption 
time constant. The time constant is a function of temperature and is given approximately by 
T = 10“ 13 exp ( Q/RT ), where Q is the adsorption energy. 
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Three equations can be written: one which relates the antechamber density of atomic oxygen to 
its sources and sinks, one relating molecular oxygen to its sources and sinks, and a similar equation for 
the wall density of atomic oxygen: 
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It is now convenient to define the time constants: 
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Assuming that dn 36 /dt is small with respect to n\^ /r 36 and dn 32 /dt is small with respect to n 32 /r 32 
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then equations (18), (19), and (20) can be rearranged and written as 
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Now, if (~A h /V) times equation (23) is added to equation (21) and if twice equation (22) is added to 
the previous sum, the following general result can be obtained: 


,16 


r i6 


16 

g 

,16 


+ 2 


32 

g 

t -32 


A w dn w 


16 


(24) 


V dt 


a g g 

The ambient atomic oxygen is the sum of a density derived from the source density of atomic oxygen, 
twice the source density of molecular oxygen, and a correction proportional to the rate of change of 
the oxygen wall density. 


Determination of the wall coverage is found in principle by integrating equation (23), given the 
four surface parameters (& 16 , a 32 , a 32 , and r^ 6 ) and values of n 36 . There are three possible sources 
for values: (.1) the mass 16 source density measurements (subtracting the portion due to the 
fractionation of molecular oxygen), (2) equation (22) (using mass 32 source density measurements 
(>ig 2 ), and (3) equation (21) (using an atmospheric model for n 36 ). The first possibility was avoided, 
because n 36 would then be the difference of large numbers and because the fractionation ratio was not 
known with complete certainty. The second possibility was initially considered the best, because it 
uses measurements having the least scatter. The third possibility suffers the disadvantage of requiring 
knowledge of n 36 , the quantity one is attempting to measure. 


An examination of the various equations shows that the surface parameters should, in principle, 
be determinable from the data. If it is assumed that adatom-gas atom recombination dominates over 
adatom-adatom recombination near perigee (as found by Wood, Baker, and Wise, 1970), derivation 
from equation (22) gives 
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and this is roughly 0.022. Also, near perigee one finds from equation (23) that & 16 is roughly equal to 
o 32 n 36 when o 32 and 1/r, 1 ., 6 are small. 
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The apogee data place further constraints on the parameters. The existence of a finite n 36 implies 
a finite t 16 , which from equation (21) is approximately 
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assuming n 36 is zero. Assuming also that the apogee mass 32 density is primarily due to adatom- 
adatom recombination, one can derive from equations (22) and (24) the values of n^ 6 at apogee, and 
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Figure 10.— Measured and calculated mass 16 
to 32 source-density ratios. 

If it is further assumed that perigee and apogee wall densities are the same order of magnitude, rough 
values of all the surface parameters will have been determined. 

An attempt was made to refine these approximate values of the surface parameters by a least 
squares fit of all the mass 16 source densities in an orbit and the mass 32 source densities at apogee. 
Values of n g 6 were predicted from equation (23), using measured n 22 and values of n ] w b calculated 
from equation (23). Values of n 22 near apogee were predicted by determining n g 6 from equation (21) 
with n 26 = 0 and substituting it into equation (22). Equation (23) was solved for n^ 6 by repeated 
numerical integration, using values of n 26 calculated from equation (22), until the n^ 6 values were the 
same at the beginning and end of an orbit. Whereas parameters could be found that improved the fit 
to the data over the rough parameters estimated above, it became apparent that a systematic variation 
in the residuals (differences between measured and predicted values) was frustrating the search for a 
satisfactory solution. 

The nature of the problem is best seen by examining the ratios of mass 16 to 32 source densities. 
A plot of these ratios versus time for a whole orbit is shown in figure 10, which is characteristic of all 
orbits after the first month. The observed ratios are remarkably symmetric about perigee and have a 
maximum at perigee. The predicted values of the ratio shown are not the best that can be achieved 
with a least squares fit, but serve to illustrate the problem. The predicted ratio has a slight dip after 
perigee and is higher before perigee than after, contrary to the near-perigee behavior of the measured 
ratio. The reason for the near-perigee behavior of the predicted ratio can be seen from equation (25). 
Assuming that molecular oxygen is the result only of adatom-gas atom recombination (independent of 
how the n ] g 6 level is established), the n g 6 /n g 2 ratio is inversely proportional to the wall density n^ 6 . 
Because the sources contributing to the wall density increase to a maximum near perigee, the wall 
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density should likewise increase toward and beyond perigee to reach a maximum where the sources 
and sinks are equal. This behavior is reflected in the predicted values of figure 10. A limiting case 
would exist if the wall density were constant, implying constant /ft 32 . In an attempt to explain the 
behavior of the oxygen ratio, a number of speculative processes were considered, including direct sur- 
face ionization of wall oxygen in the source, dissociative adsorption of 0 2 , and nonhomogeneous dis- 
tribution of ambient flux on the antechamber walls; but none was satisfactory. A process that would 
provide maximum ft 16 /^ 2 at perigee would involve the reaction of oxygen atoms (surface or gas phase) 
with some other species that peaks at perigee more strongly than atomic oxygen, resulting in a molecule 
that carries atomic oxygen into the ion source, where it is dissociatively ionized. Neither of these 
possibilities could be established from the data with certainty or used in a quantitative manner, and 
attempts to fit the data with a method that relied largely on the mass 16 to 32 source-density ratio 
were abandoned. 


The method finally used included fitting only the mass 32 source densities over the whole orbit, 
calculating / 7 32 from equation (22), using ft] 6 from equation (21) as 
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and obtaining « ^ from equation (23) as before. Lacking the measured mass 16 source density, the 
surface parameters cannot be determined without using a model atmosphere for n ] 6 . Fortunately, it 
can be shown that independent knowledge of all the surface parameters is not necessary to determine 
the density correction. Only the ratio of k l6 to a 32 is important in the correction formula if k 16 and 
a 22 are small with respect to 1 . If equation (22) is solved for ft] 6 and substituted in equation (23), 
the important terms will depend only on the k l() /a 22 ratio. Similarly, from equation (29), with 
ky 6 A /A h large with respect to 1 , the substitution into equation (23) again produces ratios of 
to a 3 ^. Thus it was decided to fix k . , at 0.022, which gives rough agreement with the perigee 

ag lo 

ft i 6 In 3 2 ratio. The value determined by Wood, Baker, and Wise (1970) was 0.044. 

g g 


It is also very difficult to separate the two wall-loss terms «J, 6 /rj, 6 - and 2a 32 (ft], 6 ) 2 because at 
apogee, where they are important, the wall density is not changing fast enough to provide a significant 
difference between linear and square-law loss. It was decided to fix 1 at 6 X 10 5 , which provides 
an approximate agreement with the apogee mass 16 source density. Numerical experiments showed 
that within wide ranges the chosen values of k l6 and r ], 6 made no practical difference to the oxygen 
correction as long as the remaining two parameters were adjusted for best fit to the mass 32 data. 

Only two surface parameters are still to be determined: o 22 and a 32 . 


The model chosen to provide the ambient atomic oxygen densities was basically the static diffu- 
sion model of Jacchia (1965), and orbits were selected that, as far as possible, had N 2 densities in 
reasonable agreement with the model. This was accomplished by comparing the exospheric tempera- 
ture, which yielded the measured N 2 density, with the exospheric temperature as predicted by Jacchia. 
The temperature used for predicting N 2 and atomic oxygen could then be corrected, to provide better 
agreement with the measured N 2 . An attempt was made, also, to pick orbits that had perigees near 
the Equator and low magnetic activity. Five orbits were selected for the determination of a 32 and a 32 . 
These occurred on Sept. 22, 24, 26, and Nov. 15, 1969, and on Mar. 22, 1970. The difference between 
the temperature needed to derive the measured nitrogen density T M and the model temperature T } is 
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Figure 1 1 .-Comparison of nitrogen ambient density with 
Jacchia model, (a) Orbit of Sept. 22, 1969. (b) Orbit of 
Sept. 24, 1969. (c) Orbit of Sept. 26, 1969. 


plotted against latitude in figure 1 1 . The approximations used for the model temperature when pre- 
dicting atomic oxygen densities are also shown. The predicted atomic oxygen densities were scaled 
with a multiplicative factor so that only the predicted oxygen latitudinal distribution and not absolute 
values would be used in determining surface parameters. This factor and the unknown surface param- 
eters were determined by a least squares fit to measured source densities for the selected orbits. 

Values of oxygen surface parameters and o^ 2 determined for the five selected orbits are given 
in table 1 . Two surface parameters determined from the orbit on Sept. 24 were chosen as the values 
to be used in ambient density processing because the N 2 densities agreed most closely with the model 
predictions, the orbit was during equinox when cross-Equator gradients are least likely, and the orbital 
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Table 1. -Surface Parameters for Oxygen a 


Date 

a 32 

aa 

32 

o° 

ag 

Sept. 22, 1969 
Sept. 24, 1969 
Sept. 26, 1969 
Nov. 15, 1969 
Mar. 22, 1970 

8.67 X 10- 19 
12.00 X 10" 19 
4.52 X 10" 19 
2.20 X 10- 19 
20.60 X 10" 19 

3.72 X 10“ 16 
3.81 X10- 16 
2.80 X 10~ 16 
1.99 X 10“ 16 
4.63 X 10“ 16 


a Based on assumed values of fc 16 - 0.022 and 1/r^f - 6 X 10 s . 


variation of densities indicated that this was an unusually quiet time period (in agreement with the 
magnetic activity indices). The comparison between the measured and predicted source densities for 
the Sept. 24 orbit is shown in figure 12; these correspond to mass 16 to 32 source-density ratios 
plotted versus time in figure 10. The predicted variation in wall density n 16 for this orbit is given in 
figure 13. 


Rather than calculate ambient densities directly from equation (29), it was decided to assume 
that was always zero at apogee and to subtract the results of equation (24) for the apogee case 
(subscript ap) from the results for other times: 
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This procedure has the virtue of insuring that the ambient density tends to zero at apogee regardless of 
small errors in calculating dn x J /dt. Perigee results will be essentially unchanged. A correction factor 
was then defined as follows: 

lit - (dn'J /*) ] 

CF = — — — (31) 

[”] 6 -(«i 6 )a P ]^ 6 + 2[«f - («| 2 )a p] /^ 32 

which means that if ambient densities are calculated by the usual method of ignoring surface effects 
but subtracting the background (the first two terms of eq. (30)), the densities corrected for surface 
effects will be (1 + CF) times the density determined from the first two terms of equation (30). A 
plot of CF versus time is shown in figure 14. 

The effect of an arbitrary change in surface parameters is demonstrated in figure 1 5 by examining 
the ratio of corrected ambient density to the model ambient density. Note the insensitivity of meas- 
ured density near 500 km on the outgoing portion of the orbit to changes in surface parameters, which 
is in agreement with the correction being very low here. 

The ambient densities derived for other orbits using the Sept. 24 surface parameters are compared 
in figure 16 with the densities derived using the best-fit parameters obtained from that particular orbit. 
There is no large discrepancy except for the Nov. 15 orbit, which, being nearer to solstice than to 
equinox, may have composition effects not accounted for in the static diffusion models. The 
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Figure 14.— Oxygen correction factor (CF) for the 
orbit of figure 12. 
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Figure 15,-Effect of surface parameter changes on 
calculated ambient density of atomic oxygen. 
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Figure 16.— Comparison of ambient densities derived from Sept. 24 surface parameters and from best fit surface 
parameters from orbits on other dates, (a) Best fit data from Sept. 22 orbit and sweep data, (b) Best fit data from 
selected orbit on Sept. 26. (c) Best fit data from selected orbit on Nov. 15. (d) Best fit data from selected orbit on 
Mar. 22, 1970. 
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comparison for the Mar. 22 equinox orbit is quite satisfactory. The fact that different pairs of values 
for and can give substantially the same ambient densities suggests that there is an effective cor- 
relation between the two parameters. This is illustrated in figure 17, where the line represents sets of 
parameters for Sept. 24 derived by arbitrarily changing one parameter and least square fitting for the 
remaining one. The parameters for the other days fall very close to this line. Thus, with the method 
used, the surface model has in effect only one degree of freedom. 

The extent to which the corrected ambient densities (determined from the measurements) depend 
on the model ambient densities used in determining the surface parameters was studied by making 
arbitrary changes in the model temperature used to predict the model density. It was found that adding 
or subtracting a constant temperature to the Jacchia model or using a constant temperature throughout 
the orbit made no significant difference in the final corrected ambient densities, although the fit to the 
source density was not always very good, and different surface parameters resulted. However, if a 
temperature gradient of ±1 K per degree latitude was added to Jacchia-model temperatures, it was 
found that this gradient was reflected in the corrected density measurements. The same happened for 
±2 K, but the degree of fit to the source densities was sufficiently poor to suggest that a gradient of 
this magnitude was unlikely in the real atmosphere for the orbit tested. The result of these tests is that 
an error in the ambient density model, which is symmetric about perigee, is unable to cause a signifi- 
cant error in the corrected measured density, but an error in the form of a temperature gradient would 
be reflected in the measured density. 

An assessment of the error resulting from the electronic failure of Aug. 13, 1969, was made by 
hand-constructing a full-orbit curve of mass 32 source density from the sweep mass 32 source density, 
as opposed to using that obtained in the step mode, and by using these data to determine ambient 
density. Presumably the sweep determinations would be more accurate at the higher altitudes. The 
densities beyond 600 km using the sweep data were systematically lower than those from the step data. 
On Sept. 22 the differences were roughly 10 percent at 600 km and 20 percent at 700 km (fig. 16(a)). 
On Sept. 26, the corresponding differences were 15 and 35 percent. In addition, the difference 
between before and after perigee densities at 600 km was about 1 0 percent smaller using the sweep 
data. 



Figure 17. -Comparison of af| and a||, assuming 
A:, 6 = 0.022 and 1/r^ 6 = 6X 10" 5 . 
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During certain time periods (e.g., Nov. 1969) there were changes near apogee in the slope of the 
source density versus time plots which suggested that some process in addition to simple desorption 
was occurring. Maximum deviations in density were on the order of 1 5 percent. Three possible 
explanations were considered: (1) temperature variation of desorption and adatom-adatom recombina- 
tion, (2) outgassing from the body or solar panels of the spacecraft, and (3) solar radiation entering 
the spectrometer antechamber. Unfortunately, all three possible effects are highly correlated in prac- 
tice, because the normal spacecraft operation was such that when the orifice was facing the Sun, the 
orifice was always facing in the direction of the solar paddles and main body. The phase was such 
that all three effects could qualitatively predict the perturbations in Nov. 1969, but would not be 
correct at other times. The ability to predict the temperature effects depends on knowledge of the 
surface interaction energies, and when these are estimated from the surface parameters already deter- 
mined, the variations in source density caused by temperature appear to be too large. Ambient densi- 
ties determined by including temperature variations differed from those determined without tempera- 
ture variations by less than 10 percent up to 600 km. There was some evidence for very small 
variations in source density during OPEP scans, which may suggest outgassing from the solar panels or 
body, but the data are too meager to give quantitative predictions for any time during an orbit. 

In summary, a set of surface parameters has been determined that provides a reasonable fit to the 
mass 16 and 32 source densities and is consistent with the ambient atomic oxygen predicted by the 
static diffusion model (Jacchia, 1965) during equinox. It was found that parameters determined from 
orbits up to 6 months apart were reasonably consistent, which indicates, in conjunction with the con- 
stant mass 16 to 32 source-density ratio at perigee, that surface conditions are relatively stable after 
the first few weeks of flight. The surface model is not applicable to the first few weeks of flight be- 
cause there was significant production of CO, and CO that has not been taken into account. 

Unfortunately the surface parameters could not be determined entirely from the source-density 
data, but depended to some degree on an assumed ambient density model. However, the gradient 
across perigee was the essential feature derived from the ambient density model, and it is probable that 
by using an orbit near equinox with perigee near the Equator the possible errors were minimized. 

During routine processing, equation (23) is solved for n by integrating numerically in serial 
fashion from orbit to orbit, without the need to repeat the process as was done for a single orbit, and 
the correction factor of equation (3 1 ) is calculated from equation (24) by use of the instantaneous 
n^ 6 . To avoid using ambient densities with large corrections, it was decided to add to the density 
error a term equal to 25 percent of the correction, if positive, and 50 percent, if negative. The dif- 
ferent error values arose because it appeared that a change in a single parameter had twice the effect on 
the correction when the correction was negative. 
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